Phase diagram of the half-filled one-dimensional t-V-V model 



O 

(N 



I 

o 



> 

(N 

od 
O 



Tapan Mishra, Juan Carrasquilla, and Marcos Rigol 
Department of Physics, Georgetown University, Washington DC, 20057 
(Dated: September 28, 2011) 

We study the phase diagram of spinless fermions with nearest- and next-nearest-neighbor interac- 
tions in one dimension utilizing the (finite-size) density-matrix renormalization group method. The 
competition between nearest- and next-nearest-neighbor interactions and nearest-neighbor hopping 
generates four phases in this model: two charge-density-wave insulators, a Luttinger-liquid phase, 
and a bond-order phase. We use finite-size scaling of the gap and various structure factors to 
determine the phase diagram. 
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I. INTRODUCTION 



The study of quantum phase transitions in strongly in- 
teracting systems has become a major field of research in 
condensed matter physics. Quantum fluctuations, rather 
than thermal fluctuations, are responsible for such zero 
temperature phase transitionsii Remarkably, quantum 
fluctuations play a more dominant role as the dimen- 
sionality of the system is reduced. In that sense, one and 
two dimensions are in general an ideal playground where 
the interplay between strong interactions and quantum 
fluctuations leads to the emergence of exotic phases and 
phenomena. For example, in one dimension, it is known 
that quantum fluctuations melt any order that would oth- 
erwise break a continuous symmetryi^i^ 

In recent years, in addition to traditional condensed 
matter systems that exhibit reduced dimensionalities, 
ultracold atomic gases in optical lattices have allowed 
the experimental realization of several strongly corre- 
lated systems that arc of great interest The degree of 
control in those experiments is such that effective one- 
dimensional regimes are now within reachi^ In addition, 
the high degree of isolation in optical lattice experiments 
enables the study of nonequilibrium phenomena not ac- 
cessible in condensed matter settingSfii^ and addressing 
questions related to the relaxation dynamics and ther- 
malization of quantum systems 

One-dimensional systems have been the focus of 
many recent works on the effects of integrability and 
phase transitions in the dynamics and thermalization of 
strongly interacting quantum systemsi^ In particular, 
one of us (in collaboration with L. Santos) studied the 
breakdown of thermalization in the t-V-V' model, and its 
relation to approaching integrable points and the various 
insulating phases present in its ground stateJ^J^ Inter- 
estingly, despite the simplicity of the model, and the fact 
that it has been studied in several previous works an 
accurate phase diagram is not available in the literature. 
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In this work, we present a detailed study of the phase 
diagram of the t-V-V model. For spinless fermions, in 
its particle-hole symmetric form, the Hamiltonian can be 



where c| and are fermionic creation and annihilation 
operators at site i, respectively. These operators obey the 
usual fermionic anticommutation relations {c^ , cj} = Si j. 

n.i = c^Cj is the number operator, t is the hopping am- 
plitude between neighboring sites, and V and V are the 
nearest- and next-nearest-neighbor interactions, respec- 
tively. 

We focus our analysis on the half-filled case, i.e., when 
the number of fermions is equal to one-half the num- 
ber of lattice sites. For V' = 0, this Hamiltonian can 
be solved exactly using the Bethe ansatzi^ As V is in- 
creased, one finds a transition from a gapless Luttinger 
liquid (LL) to a gapped charge-density wave at V/t = 2. 
For t = 0, on the other hand, the ground state of model 
(fT|) is very simple (a product state) and is determined 
by the ratio between V and V . For V' < V/2 one has a 
(... 10101010 ...) charge-density wave (CDW-I) and 
for V > F/2 one has a (. . . 1 1 1 1 . . .) charge- 
density wave (CDW-II). Here, 1 (0) denotes the presence 
(absence) of a fcrmion at a particular site. However, for 
finite values of t, a richer phase diagram emerges. 

Early numerical studies of this model were based on 
the Lanczos method^ for which only small system sizes 
can be diagonalizcd. It was found that the competition 
between t,V, and V leads to four different phases. Three 
of those phases have been already mentioned, the LL, 
CDW-I, and CDW-II phases. In addition to those, a 
bond-order (BO) phase was also found. This BO phase 
was somehow missed in some subsequent studies of this 
modelJ^ii^ More recently, the presence of the BO phase 
was discussed in a density-matrix renormalization group 
(DMRG) study, but the phase diagram was not computed 
there .-i^ 

It is interesting to note that despite the fact that the 
physics of the t-V-V model is well understood,— as men- 
tioned before, no accurate phase diagram has been re- 
ported so far. This is because large finite-size effects 
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FIG. 1. (Color online) Phase diagram of the half-filled t-V-V 
model. The phases are, from bottom to top, charge-density 
wave with a two-site unit cell, Luttinger liquid, bond-order 
phase, and charge-density wave with a four-site unit cell. The 
dashed line corresponds to V' — V/2, which, in the atomic 
limit, gives the boundary between the CDW-I and CDW-II 
phases. The hopping amplitude is set to t = 1. 



make it difficult to determine the precise boundaries be- 
tween the phases mentioned above. In this paper, uti- 
Uzing the DMRG technique, we determine the phase di- 
agram of the t-V-V model for a broad range of values 
of V/t and V /t. In order to make accurate predictions, 
we use scaling properties of several quantities, which are 
based on the universality class of the various transitions. 
Our main result is the phase diagram reported in Fig. [T] 

The exposition is organized as follows. In Sec. |TT1 we 
discuss the model under consideration and its formula- 
tion in different physical contexts. We also introduce the 
observables used to describe the different ground-state 
phases. In Sec. IIIIl we provide details on the different 
approaches used to detect the transitions from gapless to 
gapped phases. In particular, we determine the transition 
line from the LL to the CDW-I as well as the transition 
from the LL to the BO phase. In Sec. IIV| the transition 
from the BO phase to the CDW-II is studied by perform- 
ing finite-size scaling analysis of the structure factor. We 
also compare the obtained phase diagram with the one 
found in previous studies. Finally, in Sec. |Vl we briefly 
summarize our main conclusions. 



II. MODEL, OBSERVABLES, AND APPROACH 

The t-V-V model for spinless fermions [Eq. ([T])] can 
be mapped onto two other models of much interest in 
condensed matter and cold gases. 

The first of those two models is a well-known spin 
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where Sf,Sf,S^ are spin-i operators at site i, the Hamil- 
tonian (HI) takes the form 

H = Y^ [-2t (5f + SfSf^,) + VS^S^^, + VS^S!^^ 

i 

(3) 

This is the spin- 1/2 XX Z chain with an additional next- 
nearest-neighbor S^S^ interaction term. For the map- 
ping, we have assumed that the system has open bound- 
ary conditions, as will be the case throughout this work. 
For periodic boundary conditions, a boundary term ap- 
pears which depends on the total number of fermions 
(total S^) in the chain.— 

Furthermore, one can map the spin model above onto 
a model of impenetrable (hard-core) bosons, known as 
the lattice Tonks-Girardeau gas within the cold-gases 
community?^ For this, one can use the Holstein-Primakoff 
transformation for spin-I/2 particlesi^ 
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Under this transformation, Eq. can be written as 
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where al (a^ ) is the bosonic creation (annihilation) opera- 
tor obeying the bosonic commutation relations [a, , a]] = 
and n\ = a\a^ is the bosonic number operator. One 
also has an additional constraint = <\ =0 that pre- 
cludes multiple occupancies of the lattice sites. 

From the derivations above, one can see that spinless 
fermions, spins, and hard-core bosons share the same 
spectrum, and diagonal (density and 5*^) correlations. 
Hence, the phase diagram obtained in this study for spin- 
less fermions is also relevant to the corresponding spin 
and bosonic models. Experimentally, hard-core bosons 
have been realized in one-dimensional geometries in the 
presenco^i and absence^^ of a lattice. 

In order to compute the ground-state properties of the 
t-V-V model, we use the finite-size DMRG algorithm 
with open boundary conditions This method is best 
suited for (quasi-)one-dimensional problems and has been 
extensively used to study quantum spin chains.— To 
minimize finite-size effects and obtain accurate extrap- 
olations, we study systems with up to 700 sites (1000 
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in fewer eases) retaining 128 density-matrix eigenstates. 
The weight of the states discarded in the density matrix 
is less than 10~^ in all cases. To improve the conver- 
gence, at the end of each DMRG step, we use a finite-size 
sweeping procedurci^ 

To characterize the various phases of this model, we 
have studied several quantities. Here, we report results 
for the single-particle excitation gap, 

Gl = E{L, + 1) + E{L, N -I)- 2E{L, N). (6) 

It allows us to distinguish gapped and gapless phases. In 
Eq. ([6]), E{L,N) is the ground-state energy of a system 
with L sites and N fermions. 

To understand the transition to the CDW phases, we 
calculate the structure factor, which is the Fourier trans- 
form of the density-density correlation function 



Finally, the transition to the BO phase can be identi- 
fied by calculating the bond-order parameter 
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BO 



where 

In the remainder of the paper, we set t — 1. 



(8) 



(9) 



III. GAPLESS TO GAPPED PHASE 
TRANSITIONS 

We first focus on the transition between the gapless LL 
phase and the two gapped phases that surround it, which 
are the CDW-I and the BO phases. These two transitions 
can be understood using LL theory, from which one ob- 
tains that the Luttinger parameter is = 1/2 at the 
transition pointsi^ One can then use this knowledge to de- 
termine the boundaries between the LL and CDW-I/BO 
phases. The idea would be to compute K using the decay 
of correlations for finite systems, make an extrapolation 
to the thermodynamic limit, and then find the values of 
V' that for a given value of V result in K = 1/2. This 
approach was used, for example, to calculate the phase 
diagram of the extended Hubbard model.— We find that 
such a procedure leads to inconclusive results for the t-V- 
V model. The selection of the range of distances used to 
fit K from correlation functions, such as the one-particle 
density matrix pij = {c\cj), leads to a wide range of 
values of K for any given finite system. After extrapo- 
lation to the thermodynamic limit, the errors in the de- 
termination of the critical V", for each value of y, were 
found to be large. One could also try to use the appropri- 
ate functional form of the relevant correlation functions 
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FIG. 2. (Color online) Scaling of the gap LGl plotted as a 
function of V' for V = A.Q. The coalescence of different curves 
for V' > 1.0 in (a) and V' < 2.9 (b) indicates the transition 
from gapped to gapless phase. 



as obtained in low-energy effective theories, like it has 
been done in the context of spin chains2£i^ and zigzag 
laddcrsp^ but we have followed a different approach. 

Our approach is based on the study of the closing of 
the single-particle excitation gap Gl [see Eq. ©J when 
entering the LL phase. Gl can be accurately determined 
within DMRG at a relatively low computational cost (in 
particular, if one compares it with the cost of comput- 
ing correlation functions). Of course, Gl is finite for 
any finite system even in phases that arc gapless in the 
thermodynamic limit. As an example, in Fig. [5] we show 
LGl as a function of V and fixed V = 4, for two differ- 
ent intervals of V' , and various system sizes. One can see 
there that LG l is finite for all values of V' and all system 
sizes. However, in some parameter regimes LGl does not 
change with increasing system size, i.e., Gl ~ ^/L (van- 
ishes in the thermodynamic limit), while in other param- 
eter regimes, LGl increases with increasing system size 
and so Gl is finite in the thermodynamic limit. 

One can then extrapolate Gl to the thermodynamic 
limit by fitting it to a polynomial in terms of 1/L and 
obtain the value of Gl^oo by varying V' for all the values 
of V . Results for the extrapolated gap Gl^oo as a func- 
tion of V ior V = 4 are shown in Fig. [31 This plot shows 
a clear transition from a gapped to a gapless phase and 
then to a gapped phase. The system becomes gapless at 
V^' 1.0 and the gap reopens at V ^ 2.9. Still, within 
this approach, one encounters the difficulty of pinpoint- 
ing the exact values of V' for the phase transition. 

In order to obtain accurate values of V for the bound- 
aries between the gapped and gapless phases, we make 
use of the knowledge that the transition between the 
gapless and gapped phases belongs to the Berezinskii- 
Kosterlitz-Thouless (BKT) type,^ At the BKT transition 
the gap closes as 
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where a is a constant. 

The correlation length which is closely related to the 
gap, is finite in the gapped phase, diverges at the critical 
point as ^ G~^, and remains infinite in the LL phase. 
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FIG. 3. The extrapolated gap in the thermodynamic limit is 
plotted as a function of V' for V — 4. There is a transition 
from a gapped to a gapless phase at V' ~ 1.0 and then to a 
gapped phase again at V' ~ 2.9. 



Wc utilize the following finite-size-scaling relation for the 
gap in the vicinity of the phase transition, 



LGlxIi 



1 



2 In L + C 



F 



(11) 



where is a scaling function and C is an unknown 
constant to be determined. This scaling ansatz resem- 
bles the analogous scaling relation for the resistance in 
the charge-unbinding transition of the two-dimensional 
Coulomb gas.— In such a transition, which is also of the 
BKT type, the resistance near the critical point behaves 
as the gap does in the gapped to LL transition; namely, 
it vanishes exponentially as in Eq. pO)) . 

At the critical point and close to it within the LL phase, 
one expects the values of F{^/L) to be system-size in- 
dependent because of the divergence of the correlation 
length; i.e., plots of LG'^ = LGl [1 + l/{2lnL + C)] as 
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FIG. 4. (Color online) The scaled gap LG'j^ as function of the 
XL ~ InL — a/ \/\V' — V^\- The insets show the scaled gap 
as function of V' . (a) Transition from the CDW-I to the LL 
phase, (b) Transition from the LL to the BO phase. Both 
transitions are presented for a fixed value oi V = 4. 



function of V for different system sizes should merge in 
that region. Furthermore, the resulting curves of plotting 
LG'j^ as a function of ^/L for several values of L should 
also be system-size independent. Equivalently, one can 
plot LG'j^ as a function of xl = InL — In ^ to get the 
collapse of the curves. In our calculations, the values 
of a, C, and arc fitted to produce the best possible 
collapse of the data within the gapped side of the phase 
transition. On that side, wc take into account that the 
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correlation length diverges as ^ ^ exp 

We have tested the accuracy of this procedure by deter- 
mining the critical value of V for the transition between 
the LL and CDW-I phases when V = 0. As mentioned 
in the introduction, the critical value of V in this case is 
known analytically (the Heisenberg point). We consid- 
ered systems with up to 700 sites and obtained a critical 
value Vc = 2.02 ± 0.01, which is consistent with the ana- 
lytical result. 

When V' 0, this model is not exactly solvable. In 
Fig. m we show results for LG'j^ vs xl (main panels), as 
well as LG'j^ vs V (insets), for y = 4 and several val- 
ues of L. The main panels make evident the collapse 
of all the data for the rcscalcd gap as a function oi xl, 
shown here only within the gapped phases. Two impor- 
tant limits can be understood by analyzing the collapse 
curves. As the critical point is approached {xl ^oo), 
the scaling function approaches a constant value, which 
in turn, implies the vanishing of the gap. In the limit of 
large x^, one can see that the scaling function increases 
rapidly, which is necessary for the gap to be finite in the 
gapped phases. We then confirm two aspects of the BKT 
transition, the exponentially divergent correlation length 
and the logarithmic corrections to the gap. The insets in 
Fig. |4] clearly show that, when coming from the gapped 
phases, the LG'j^ curves merge at a point [V^ = 1.16 for 
the CDW-I to LL phase, Fig.gja), and = 2.55 for the 
BO to LL phase, Fig.m^b)] and remain close to each other 
for some finite region within the gapless phase. Moreover, 
as the size of the system increases, the interval over which 
the curves are seen to collapse within the LL phase in- 
creases. As said before, because of the divergence of the 
correlation length in the LL phase, the curves should be 
system-size independent in the vicinity of the transition 
point, where the scaling relation ([TT]) holds. 



A. CDW-I to LL transition 

As mentioned in the introduction, only two phases are 
present in the system for V' ~ 0. Those are a LL phase 
for < < 2 and a CDW-I phase for V > 2. By adding 
a small V , the transition between LL to the CDW-I is 
shifted to larger values of V, thus enhancing the stability 
of the LL phase. This is because the interplay between 
CDW configurations of the type (... 101010 ...), 
favored by increasing values of V, competes with those 
of the type (... 11001100 ...), favored by increasing 
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FIG. 5. (Color online) Density-density structure factor 5(71) 
as a function of 1/L for different values of V' , for V = 4. 
S{n) is predicted to be finite in the thermodynamic limit for 
V' < 1.2, making evident the existence of a CDW-I phase in 
that region of parameter space. 

values of V' , thus making the system gain kinetic energy 
because of the net reduced charge ordering effects. Nev- 
ertheless, for large enough values of V, the CDW-I phase 
is prevalent. In order to see that this is the case, we 
have calculated the density-density structure factor S{k) 
as defined by Eq. ([7]). Because of the charge order, the 
CDW-I phase is characterized by a finite value of S{tt) 
in the thermodynamic limit. In Fig. [Sj we show S{7r) as 
function of 1/L for different values of V , and for V = 4. 
It is apparent in that figure that, for V' < 1.0, S'(7r) 
is finite in the thermodynamic limit. For V' > 1.2, on 
the other hand, S{tt) becomes very small as the system 
size increases, which suggests that it will be zero in the 
thermodynamic limit. This is consistent with the results 
from the scaling of the gap that, for V = 4, predicted 
the critical value of V for the transition between the LL 
and the CDW-I phase to be V' = 1.16. Despite the fact 
that with the finite-size scaling of the structure factor 
S{t:) it is difficult to locate the exact transition point, 
this calculation provides evidence of the nature of the 
CDW-I phase and the way the charge order develops as 
one crosses the critical region. 



FIG. 6. (Color online) Finite-size scaling of Obo for different 
values of V', for V — 4, across the LL to BO phase. 

by the Obo parameter, the system exhibits no charge 
order. Therefore, two quantities can be used to charac- 
terize the BO phase, the gap and Obo, which are both 
finite in that phase. In a finite-size system, however, BO 
oscillations are present in the LL phase, but vanish as 
the system size in increased. 

Following the arguments in Ref.lsolandlsil the strength 
of such oscillations is assumed to decay as . This 
scaling relation holds inside the LL region and in partic- 
ular at the transition point between the LL and the BO 
phase, where the Luttinger-liquid parameter takes the 
value K = 1/2. Because of that, we extrapolate the BO 
parameter using L~^/^ rather than L~^. In Fig. [6l we 
plot Obo as a function of L~^/^ for different values of V' 
a,t V = 4. One can see there that Obo extrapolates to a 
nonzero value in the thermodynamic limit for V' > 2.6. 
Despite the large size effects that are present at the BKT 
transition, for values of y < 2.5, Obo extrapolates to 
very small values. These results are compatible with the 
critical value ~ 2.55 obtained from the finite-size- 
scaling analysis of the opening of the charge gap. Note 
that in Fig. El as the values of V decrease and approach 
the critical value 14, the curves become closer to straight 
lines. This supports our assumption that, at the transi- 
tion point, the BO oscillations decay as L~^/^. 



B. LL to BO transition 

At the other boundary of the LL phase, when V' > 
V/2, there is an instability toward the formation of 
a gapped BO phaseJ^^ In the limit of vanishing V, 
the bond-order phase arises because of the competition 
between the kinetic energy term and the next-nearcst- 
neighbor interaction that induces a CDW-II order. A fi- 
nite value of V competes with the bond ordering induced 
by V' , thus enhancing the stability of the LL phase and 
increasing the critical value of V. As mentioned pre- 
viously, the BO phase is characterized by a finite value 
of Obo [see Eq. ([5|)] in the thermodynamic limit. Al- 
though this phase exhibits bond oscillations as unveiled 



IV. GAPPED TO GAPPED PHASE 
TRANSITION 

A. BO to CDW-II transition 

Ultimately, when V' is very large, the system always 
forms a CDW-II insulator as the energy of the configura- 
tion (... 110011...) becomes energetically more fa- 
vorable. Hence, within the gapped region with V' > V/2, 
there is an additional phase transition between the BO 
and CDW-II phases. The latter phase is characterized 
by a finite value of S{'k/2) in the thermodynamic limit. 

Taking S{'k/2) as the order parameter for the CDW-II 
phase, we obtain the critical point accurately by means 
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FIG. 7. (Color online) Scaled S{n/2) vs the scaled control 
parameter for V — 4. The collapse of all the data points onto 
a single curve confirms the critical point at — 4.293. Inset: 
Crossing of the scaled S{n/2) for V = 4 and different system 
sizes. All the curves intersect at = 4.293 indicating the 
transition point. 



of scaling theory. We start with the ansatz 

5(7r/2)L2^/'' = F ((1/' - ]/') L^/") (12) 

where /? is the critical exponent corresponding to the uni- 
versality class of the phase transition and v is the correla- 
tion exponent, is a scaling function and is now the 
critical point for transition between the BO to CDW-II 
phases. 

Comparing results for different system sizes as V' is 
changed, we found that the curves for different system 
sizes intersect at a unique critical point if /3 = 0.125 and 
ly — 1 (these are the exponents of the 2D Ising univer- 
sality class). Results for a value of 1^ = 4 are shown 
in Fig. El for which ~ 4.293. The scaled S{tt/2) as 
a function of scaled control parameter V' is plotted in 
Fig. [71 The collapse of the data confirms the accuracy of 
as determined from the crossing in the inset of Fig. [71 
Notice that, as V' is increased, the transition between 
the BO to CDW-II occurs while the charge gap contin- 
ues to grow monotonically, i.e., it does not vanish at the 
transition point. See the behavior of the charge gap in 
Fig. |3lfor values around the critical ~ 4.293. 

In order to have an independent check that the previ- 
ous approach leads to the correct results, we have also ex- 
trapolated the value of S{tt/2) to thermodynamic limit. 
The results from this procedure, once again for V ^ A, 
are presented in Fig. El There, we plot S{tt/2) vs l/L 
for different values of V' . The values of the extrapolated 
S{tt/2) clearly converge to zero for V < V^. 

A complete analysis like the one presented so far, but 
for different values of V, allowed us to determine the 
phase diagram presented in Fig. [H 

Finally, to allow for a direct comparison with early 
calculations for this model, reported in Refs. [l^ and [13, 
we present in Fig. [HI the phase diagram in the plane 
t/V — V' /V. This phase diagram is in qualitative agree- 



FIG. 8. (Color online) Scaling of S{n/2) vs l/L for F = 4 
and different values of V' between V' = 4.25 and V' = 4.32 
in steps of 0.01. One can see that S{tt /2) extrapolates to zero 
below the critical value Vl = 4.293. 
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FIG. 9. (Color online) The same phase diagram as in Fig. [T] 
but now in the t/V vs V' /V plane. 



mcnt with the one reported in Refs. |l2| and |l7|. How- 
ever, quantitative differences on the precise location of 
the phase boundaries are apparent. Although we did not 
study the region of the phase diagram where t/V — 5- 0, 
the trend of our transition lines suggests that they will 
cross for small but finite values of t/V ^ implying a tri- 
critical point J^iii 



V. CONCLUSION 

We have presented a comprehensive study of the phase 
diagram of the t-V-V model in one dimension. Using the 
scaling of the gap in the BKT transition between the LL 
phase and the CDW-I/BO phases, we have obtained ac- 
curate results for the boundaries between the gapless and 
gapped phases of this model, which are confirmed by the 
extrapolation of various order parameters. The phase 
transition between the BO phase and the CDW-II phase 
(both of which are gapped) was also studied using scaling 
theory for the density-density structure factor. This lat- 
ter phase transition was found to be much less sensitive 
to finite-size effects and the exponents computed were 
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